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Transiting exoplanets from the CoRoT space mission 

V. CoRoT-Exo-4b: Stellar and planetary parameters * 
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ABSTRACT 

Aims. The CoRoT satellite has announced its fourth transiting planet (Aigrain et al. 2008) with space photometry. We describe and analyse 
complementary observations of this system performed to establish the planetary nature of the transiting body and to estimate the fundamental 
parameters of the planet and its parent star. 

Methods. We have analysed high precision radial- velocity data, ground-based photometry, and high signal-to-noise ratio spectroscopy. 
Results. The parent star CoRoT-Exo-4 (2MASS 06484671-0040219) is a late F-type star of mass of 1.16 M G and radius of 1.17 R B . The planet 
has a circular orbit with a period of 9.20205d. The planet radius is 1.19 Rj up and the mass is 0.72 Mj llp . It is a gas-giant planet with a "normal" 
internal structure of mainly H and He. CoRoT-Exo-4b has the second longest period of the known transiting planets. It is an important discovery 
since it occupies an empty area in the mass-period diagram of transiting exoplanets. 
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1. Introduction 

The masses and radii of extrasolar planets are fundamental pa- 
rameters for constraining the physics at play in their interiors. A 
combination of photometric and Doppler techniques is presently 
the only way to indirectly derive these parameters. However, this 
possibility is limited to planets that transit in front of their stars, 
hence a sample of objects that is highly biased towards close-in 
orbits. 

The planet search programme conducted from space by 
CoRoT (Baglin et al. 2006) is currently monitoring tens of thou- 
sands of stars and is regularly discovering new transiting plan - 
ets dBarge et al.ll2008l: lAlonso et al.ll2008l: IPeleuil et alj |2008a). 
CoRoT has the ability to discover planets with smaller radii 
and/or longer periods than what is routinely achieved by ground- 
based transit surveys. The evolution of fundamental parameters 
of extrasolar planets will thus be explored as a function of the 
distance to the star and as related characteristics like irradiation. 
This has important implications for our understanding of these 
planets. For example, it is a way to test whether tidal damp- 
ing (highly dependent on orbital distance), increased opacities, 
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or alternative theories are the explanation f or the anomalously 
large sizes of some extrasolar planets (e.g. iGuillotetalJ 12006: 
iBurrows et ai]|2007t IChabrier & BarafFe 2003). 



With the current detection threshold of the first CoRoT 
run (45 day duration), 26 transiting events were detected. Two 
of them tu r ned out to be plan e ts ( CoRoT-Exo-lb and 4b, 
iBarge et alj d2008l) : I Aigrain et all d2008l) . 2 are grazing eclips- 
ing binaries, 7 are low-mass binaries, and 10 are background 
eclipsing binaries. Four events remain unsolved cases, either be- 
cause they are fast rotators (no possible radial velocity follow- 
up, the cross-correlation function is too wide or undetected with 
HARPS) or because they showed no radial- velocity variations. 
These cases are still in the list of systems to be observed at 
the next visibility season with higher accuracy. The statistics 
achieved on this run is compatible with expectations and the 
simulations on CoRoT detection, with its current threshold and 
follow-up performances. 

The present paper focuses on the fourth planet discovered by 
CoRoT. The detect ion and light curve analysis is described in a 
companion paper ( Aigrain et al.l2008l hereafter Paper IV). With 
its 9.20205d period, the planet CoRoT-Exo-4b has the second 
longest period kno wn so far in the fam ily of transiting planets, 
after HD 17156b dBarbieri et al.ll2007l) . We conducted a series 
of complementary observations from the ground to establish the 
planetary nature of CoRoT-Exo-4b and to derive the fundamen- 
tal parameters of both the star and the planet. 
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2. Observations 

2.1. Photometric observations 

The photometric measurement was performed on-board CoRoT 
in an aperture that is adapted to each star and is trapeze shaped, 
with a solid angle close to that of a 14 arcsec diameter cir- 
cle. It means that a faint eclipsing binary in the vicinity of the 
main target can produce a transiting event diluted enough to 
mimic a planetary signal dBrownl l2003). The contamination fac- 
tor of CoRoT-Exo -4 is estim ated from the CoRoT entry cata- 
logue EXODAT (lDeleuil et al.ll2008bl) to be extremely low, i.e. 
only 0.3% of the light comes from neighbouring stars. The first 
step in establishing the nature of the transiting body is to check 
the photometric behaviour of the stars in the vicinity of the tar- 
get, and especially those that fall inside the photometric aperture 
mask. This is easily done with ground-based telescopes. 

CoRoT-Exo-4 was first observed fr om the Wise Observatory 
18-inch telescope dBroschet al.ll2008h on November 12, 2007. 
Only a partial transit was recorded, thus preventing a precise es- 
timate of the transit depth (of about 1 .4%, which should be com- 
pared to the out-of-transit scatter of 0. 16% in 5-min bins), but it 
is important that no significant variation was observed in nearby 
stars. Another transit event was observed with the Euler Camera 
(120-cm Swiss Telescope at La Silla Observatory) on November 
30th, 2007. Two time series of about 40-min were collected in 
the star field, with a reference series taken at a similar elevation 
but one night prior to the transit. The transit depth was estimated 
to be 1.2 + 0.15%, in agreement with the CoRoT measurements 
(1.1%). Finally, the ingress of another transit was observed at the 
Canada-France-Hawaii Telescope with MegaCam on December 
10, 2007. Again, only a partial transit was observed, and the tran- 
sit depth was not fully reliable. 

The three sets of photometric follow up confirm that 1) the 
transit observed by CoRoT occurs on the central target at the 
predicted epochs, 2) no star encompasses a deep eclipse, or even 
a short-term variation, in the vicinity of the target, 3) the transit 
depth is compatible with the depth measured by CoRoT (i.e. no 
strong dilution by background objects). A detailed analysis of 
these data will be given in Deeg et al. (in preparation). 



2.2. Radial velocity observations 

Seven radial velocity measurements of CoRoT-Exo-4 were ob- 
tained with SOPHIE on the 1 .93m telescope at Observatoire de 
Haute Provence (France), from October 31, 2007 to February 2, 
2008. These first data points established the planetary nature of 
the companion, showing a detectable and low-amplitude, radial- 
velocity variation in phase with the CoRoT ephemeris. Thirteen 
additional measurements were also obtained with HARPS (3.6m 
telescope at La Silla, ESO, Chile), from November 19, 2007 to 
January 21, 2008. 

Both instruments are cross-dispersed, fiber-fed echelle spec- 
trographs on which precise radial-velocity measurements were 
performed with regular Th-Ar calibration observ ations and the 
reduction involves c ross-correlation techniques (Baranne et al.l 
ll996l:lPepeetal.l2002l) . With SOPHIE, the high-efficiency mode 
was used with a resolving power of 40,000, and the sky back- 
grou nd was simulta neously recorded to correct for contamina- 
tion dBouchvl 12006). The radial velocities were obtained with a 
weighted cross correlation of the extracted spectra with a numer- 
ical mask corresponding to a G2 star. 

The data are listed in Table [2] (available in the on-line ver- 
sion) and shown in Fig. Q] The orbital solution was derived us- 
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Fig. 1. Radial-velocity measurements of CoRoT-Exo-4 as a 
function of time, showing the best-fit solution and O - C residu- 
als. Black filled circles are HARPS data, and red open diamonds 
are SOPHIE data. 



ing the combined HARPS and SOPHIE datasets. In the fit we 
included a constant offset between the two sets of measure- 
ments. The ephemeris of the CoRoT light curve was used to 
constrain the solution. The data were best-fitted with a circular 
Keplerian orbit of semi-amplitude of 63 +6 m/s. When using the 
13 HARPS data points alone, the residual rms is 17.6 m/s. With 
the combined SOPHIE and HARPS datasets (20 measurements), 
the residual rms is 19.4 m/s. The offset between SOPHIE and 
HARPS data is only 5 m/s. Figures Q] and [2] show the data and 
the best fit as a function of time and orbital phase. The O - C 
time series does not show any variations that would indicate a 
second massive body in the system with a baseline of 100 days 
(Fig.[D. 

The bisector analysis is shown in Fig. [3] The errors on the 
span of the bisector slopes are twice as large as for the peak 
position of the cross-correlation function. The behaviour of the 
bisector span (BIS) with respect to the radial velocity shows a 
dispersed relation with little or no trend. Given that short-term 
variations are seen in the CoRoT light curve (see with Fig. 2 in 
Paper IV), we tentatively interpret these variations as possibly 
stemming from a residual of stellar activity. 

2.3. Determination of spectroscopic parameters 

From radial velocity observations we estimated the projected ve- 
locity of the parent star. The result is vsin/ = 6.4 + 1 km/s, 
based on the width of the cross-correlation functions from 
SOPHIE and HARPS. This agrees with a direct comparison 
of the observed spectrum with a rotationally broadened syn- 
thetic spectrum. With our determination of the stellar radius 
Ri, = 1.157? (see below), this translates into a rotational pe- 
riod of the star ranging between P rot = 8.5 and 10.0 days. This 
value for the rotation period is compatible with the analysis in 
Paper IV, and is consistent with a synchronicity of the stellar 
rotation and the planet orbit. 
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Fig. 2. Radial-velocity measurements of CoRoT-Exo-4 with 
respect to the orbital phase as estimated from the CoRoT 
ephemeris and the best-fit solution. Black filled circles are 
HARPS data, and red open diamonds are SOPHIE data. 
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Fig. 3. Bisector variations (span of the bisector slope) as a func- 
tion of orbital phase (top panel) and radial velocity value (bot- 
tom panel). Black filled circles are HARPS data, and red open 
diamonds are SOPHIE data. 



We determined the fundamental parameters of the host star 
by analysing a HARPS spectrum. The signal-to-noise (S/N) level 
in the continuum is around 60 per pixel in the range 5 000- 
6 000 A and it decreases to 40 in the red and blue parts. The 
spectrum was carefully normalized by identifying continuum 
windows in a synthetic spectrum calculated with the parameters 
of the target star. A spline function was fitted to the continuum 
windows and divided into the observed spectrum to obtain the 
normalized spectrum. We made sure that the shapes and depths 
of lines in adjacent spectral orders were in good agreement. In 
the final spectrum, the orders were merged and the overlapping 
parts added using weights based on the measured S/N. 



We analysed the normalized sp ectrum using the semi - 
automatic software package VWA (IBruntt et all 120041 120081) . 
More than 600 mostly non-blended lines were selected for analy- 
sis in the wavelength range 3 980-6 830 A. The VWA uses a tmo- 
sphere models interpolated in the grid by Heiter et al. (2002) and 
atomic parameters from the VALD database (Kupka et al. 1999). 
Abundances were c alculated relative to the solar spectrum from 
Hink le et al] (J2000) following the approach by Bruntt et al] 
(2008). The VWA automatically adjusts the parameters (T e g, 
log g, microturbulence) of the applied atmospheric model to get 
consistent results for all iron lines. The program iteratively min- 
imizes the correlation between the abundance and equivalent 
width found from weak Fe i lines (equivalent widths < 80 mA), 
the correlation of abundance and excitation potential of Fe i lines 
(equivalent widths < 140 mA), and the mean abundance of neu- 
tral and ionized Fe lines must yield the same result. To estimate 
the uncertainty on the derived parameters and abundances of 
individual elements, we measured the sensitivity of the results 
when changing the parameters of the atmospheric model (see 
IBruntt et aljl2008l) . 

The result is r eff = 6190+60 K, logg = 4.41 ±0.05, and a mi- 
croturbulence £, = 0.94 + 0.07 km s~'. The determined spectro- 
scopic parameters do not depend on the adopted value of v sin i. 
The overall metallicity is found as the mean abundance of the 
elements with at least 20 lines (Si, Ca, Ti, Cr, Mn, Fe, Ni) yield- 
ing [M/H] — +0.05 + 0.07. Details on individual abundances of 
CoRoT host stars will be given in a later paper. 



3. Estimation of the stellar and planetary 
parameters 

3.1. Stellar mass, radius, and age 



Using evolutionary tracks (ISiess et al.ll2000t iMorel & Lebretor] 
2007), the spectroscopic determination of T e ff and the ratio 

1 /^ 

M s IRs estimated from the CoRoT photometry (Paper IV), we 
derived the following estimates for the mass and radius of the 
star: M, = 1.16*.9;9JM andfl, = 1.17*9/9J.R . This corresponds 
to a surface gravity of log g - 4.37 ± 0.02, which agrees with the 
spectroscopic value of 4.41 + 0.05. 

While the details of the process of how lithium is depleted in 
F-stars is not yet fully understood, it is now well-established that 
lithium is depleted for stars w ith a temperature below 6 200 K 
(Steinhauer& Deliyannis 2004). However, lithium depletion in 
F-stars is expected to be slow, and data taken for stars with the 
same age and spectral type show a huge scatter in the equiva- 
lent width. By comparing the equivalent width of the Li i line 
at A6708 (3 5.2 ± 4.8 mA) of CoRoT-Exp-4b with the data fo r 
the Hyades (Boesgaard & Tripicco 1986; Thorburn et al. 1993), 
the UMa moving gro up (Soderblom et alll993l) . and NGC 3960 



(Prisinzano & Randich 2007), an age of about 1*1-9 Gyr is the 



most likely estimate for CoRoT-Exo-4b. 



3.2. Planetary mass and radius 

Using the values obtained from the transit fit to the CoRoT light 
curve (Paper IV), the characteristics of the radial-velocity fit, 
and the stellar parameters from the previous section, we get 
a planetary mass of M p - 0.72 + 0.08 Mj up and a planetary 
radius R p - 1.1 9*9-95 /?j U p, The semi-major axis of the orbit 
is a = 0.090 + 0.001 AU. The mean density of the planet is 
0.525+0.15 gem 3 . This new plan et is found in the ex pected lo- 
cation of the mass-radius diagram (Fre ssin et alj|2007l) . In terms 
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Table 1. The orbit parameters and the fundamental parameters 
of CoRoT-Exo-4b and its host star. 



age (Ga) 



Fig. 4. Theoretical models for the contraction of CoRoT-Exo- 
4b with time com pared to the observations. The models 
JGuillot et alj 2006) assume either solar composition (left panel) 
or the presence of a central dense core of 20 M e (right panel). 
In each case, three possible evolution models are shown: (1) a 
standard model; (2) a model assuming that 1 % of the incoming 
stellar luminosity is transformed into kinetic energy and dissi- 
pated in the deep interior by tides, (3) a model in which interior 
opacities have been arbitrarily increased by a factor 30. 



of composition, Fig. [4] shows that it contains a relatively small 
amount of heavy elements (0 to up to ~ 30 M ffi , depending on 
the model chosen), as expected from the empirical correlation 
betwee n stellar metallicity and mass of heavy elements in th e 
planet dGuillot et alJl2006t iBurrowset al]l2007t lGuillotll2007h . 
The agreement with different models of planetary evolution im- 
plies that we cannot rule out any specific scenario used to ex- 
plain anomalously large planets at this point. However, the dis- 
covery of other transiting planets on long-period orbits will be 
extremely valuable. 



4. Discussion 

The planet CoRoT-Exo-4b has the second longest period of the 
transiting systems known today and is found in a region of the 
mass-period parameter space that was previously empty (see 
Fig. [5] available on-line). 

The mean density of 0.525 g cirT 3 makes CoRoT-Exo-4b a 
gas-giant with a "normal" internal structure dominated by H and 
He. It is located in the bulk of Hot Jupiters in the mass-radius 
diagram. We estimate the age of the parent star to be about 1 Gyr; 
thus, the star is likely to be fairly active and its age and rotation 
rate (Paper IV) are compatible. 

New discoveries made with CoRoT will likely populate the 
mass-period diagram in the period range 5-50 days and allow 
further investigations of the importance of star-to-planet distance 
on the internal structure of giant planets. This will complement 
the output from current radial-velocity surveys, allow us to better 
understand the mechanisms that govern the contraction of these 
planets, and provide constraints for planet formation models. 
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32.301 ± 0.005 
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Table 2. Radial-velocity measurements of CoRoT-Exo-4 with SOPHIE (lines 1 to 7) and with HARPS (lines 8 to 20). "BIS" is the 
slope of the bisector span, measured in the cross-correlation function. 



BJD 

- 2400000 



RV 

[kms- 1 ] 



Uncertainty 
[kms- 1 ] 



BIS 

[km s -1 ] 



err(BIS) 
[kms- 1 ] 



54405.6545 


32.370 


0.024 


-0.021 


0.048 


54408.6667 


32.278 


0.018 


-0.004 


0.036 


54410.6344 


32.258 


0.027 


-0.025 


0.054 


54413.6784 


32.343 


0.034 


-0.113 


0.068 


54416.6908 


32.294 


0.027 


-0.097 


0.054 


54494.3741 


32.259 


0.019 


+0.020 


0.036 


54507.3778 


32.320 


0.020 


+0.067 


0.040 


54423.8430 


32.352 
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